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ABSTRACT 

Within the next few years gravitational waves (GWs) from merging black holes (BHs) and neutron 
stars (NSs) may be directly detected, making a thorough theoretical understanding of these systems 
a high priority. As an additional motivation, these systems may represent a subset of short-duration 
gamma-ray burst (sGRB) progenitors. BH-NS mergers are expected to result from primordial, quasi- 
circular inspiral as well as dynamically formed capture binaries. The latter channel allows mergers 
with high eccentricity, resulting in a richer variety of outcomes. We perform general relativist ic 
simulations of BH-NS interactions with a range of impact parameters, and find significant variation 
in the properties of these events that have potentially observable consequences, namely the GW 
signature, remnant accretion disk mass, and amount of unbound material. 

Subject headings: gravitation — gravitational waves — black hole physics — stars: neutron — Gamma-ray 
burst: general 



1. INTRODUCTION 

Merging binaries consisting of black holes (BHs) and 
neutron stars (NSs) are prime targets for observation by 
ground-based gravitational wave (GW) detectors (such 
as LIGO, Abramovici et al.llT992.) and may be the pro- 
genitors of some sho rt-hard gamma-ray bursts (sGRBs, 
iNaravan et al. I [19921 ). The great diversity of sGRB char- 
acteristics and the potential variation in the correspond- 
ing GW signals motivates a thorough investigation of 
the possible outcomes of binary compact object (BCO) 
mergers. BCOs may form through evolution of primor- 
dial binaries or through dynamical process es in star clus- 
ters (jO^Learv et al.l I2QQ9I : iLee et al.l lMlQ). The latter 
population motivates this study of BH-NS interactions 
for systems which are initially marginally unbound. 

Star clusters at the centers of galaxies undergo mass 
segregation, resulting in heavier objects concent rated to- 
ward the center (see, e.g. iBahcall fc W^Illl977[ ). Recent 
Fokker-Plank models suggest that the Galactic nuclear 
cluster (NC) should have ^ 1800 BHs and ^ 4 00 NSs 
in the central 0.1 pc (i Hopman fc A lexander '20061). Such 
clusters are thus promising sites for BH-NS close encoun - 
ters. Using models of galactic NCs. lO'Learv et aTl (|200^) 
calculate the rate of binary BH formation through GW 
emission in close encounters. They find corresponding 
Advanced LIGO detection rates between 5 and 2700 per 
year, and estimate that the BH-NS rate could be around 
1% of this. These capture binaries form with relatively 
small periapsis separations, r^; in particular, ~ 30% form 
with Vp < lO M, where M is the to tal mass of the sys- 
tem (Fig. 4 of lO'Learv et al.l (j2009f )). (Unless otherwise 
stated we employ geometric units with G = c = 1.) This 
is due in part to the large velocity dispersion in the clus- 
ter core (^ 1000 km s~^), but also to gravitational fo- 
cusing, which may be understood as follows. The to- 
tal rate is proportional to the cross-section: F oc 7r6^, 
where h is the impact parameter. However, for Newto- 



nian hyperbolic orbits with relative velocity w at infinity, 
Vp = \?'uP' I2M + 0(w^\ Thus, the rate is linearly pro- 
portional to Tp. 

Globular clusters (GCs) that have undergone core col- 
lapse may also host BCO close encounters d ue to the high 
density of compact objects in their cores (j Fabian et al.l 
119751 : iGrind lav et al. 2006). For example, models of M15 
calibrated to the observational velocity dispe rsion yield 
a NS fr action of ^ 55% in the inner 0.2 pc (|Dull et al.l 
Il997[ ). iLee et all (|2010f ) calculate the expected rate of 
BCO interactions inside Ml 5 as a function of time and 
then scale these results for GCs with a distribution of 
half-mass relaxation times. Depending upon the GC evo- 
lution model, they find that the global rate for BH-NS 
collisions (i.e., events for which < R^s + ^bh) peaks 
at 8 — 25 yr~^Gpc~^ at redshifts between z = 0.36 
and z = 0.97, and slowly declines to between 50 — 85% of 
peak by z = 0. (We obtained the BH-NS collision rate by 
re-scaling their NS-NS res ults according to the factors in 
Table 3 of ILee et al.ll2010D. F or the fiducial BH-NS sys- 
tem considered bv iLee et al.l (|2010f ). collisions occur at 
Tp < 2.7 M. Since the rate scales linearly with r^, this 
implies an interaction rate ~ 30 — 100 yr~^Gpc~^ with 
(for example) Vp < 10 M. 

Population synthesis models (jBelczvnski et al.1 l2010f ) 
find comparable rates for primordial BH-NS merg- 
ers: from ^ 0.1 yr ~^Gpc~^ (pessimistic) to ~ 
120 yr ~^Gpc~^ (optimistic). However, primordial BH- 
NS binaries will e nter the LIGO band w ith essentially 
zero eccentricity (|Kowalska et al.l l2011f ). Thus, GW 
signals from BH-NS close encounters should be read- 
ily distinguishable due to their significant eccentricities. 
We further note that the sGRB rate, 8-30 yr~^Gpc~^ 
(|Guetta fc PiranI [M06), is comparable to the primor- 
dial BH-NS rate, and somewhat less t han the expected 
NS-N S merger rate (30-40 yr~^Gpc~^ (iBelczvnski et~all 
I2010f )). The estimates of ILee et al.l (|2010D thus suggest 
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that close encounters in clusters could contribute signif- 
icantly to the sGRB progenitor population, especially if 
their emissio n is less tightly beame d than that of primor- 
dial mergers (Grindl av etldl [20061 ). 

BH-NS mergers in clusters with < lOM will exhibit 
complicated behaviors probing the strong field regime 
of general relativity (GR). For stellar mass BH com- 
panions, one cannot treat the NS as a perturbation of 
the BH spacetime. Furthermore, the non-linear nature 
of GR will most strongly manifest during a close en- 
counter of the BH-NS. Numerical simulations within full 
GR are thus the preferred tool for exploring these sys- 
tems. To date, such si mulatio ns have been performed by 
several grou ps (Etien ne et al.l[2 009: Kvutoku et al. 2010; 
Duez et al] p 008: Ch awla et al.i i2010» : iPannarale et al.l 
201 1[ ) fseelDuez 2010 for a review). These studies employ 
quasi-circular initial data, appropriate for primordial sys- 
tems, and have explored a range of behaviors depending 
on the mass ratio, the BH spin, the NS equation of state 
(EOS), and the magnetic field. Our results complement 
these works by offering a first study within full GR of ini- 
tially hyperbolic encounters, of relevance to BH-NS cap- 
ture events and related systems that merge with large 
eccentricity. 

In the remainder of the letter we outline our numerical 
method, present results of our parameter space survey, 
and discuss some of their implications. The main result is 
the striking dependence of the outcome — disk mass, un- 
bound material, and GW signal — on the impact parame- 
ter. Though the most "extreme" outcomes might require 
fine-tuning and hence be rare, there is strong variation 
over the entire range Vp < lOM we considered: for a hy- 
perbolic encounter in a NC, this corresponds to roughly 
30% of encounters that lead to a bound system. We were 
not able to follow the larger cases through merger due 
to lack of computational resources, though certainly a 
fraction of these should also exhibit similar variability 
at the time of merger. This suggests that these systems 
could be a wellspring of varied and interesting GW and 
electromagnetic emission. 

2. NUMERICAL APPROACH 

Our 3D numerical code solves the Einstein field 
equations using finite difference (FD) techniques with 
Berger and Oliver style ad aptive mesh refinement 
(AMR, iBerger fc Oligerl Il984f ). The numerical scheme 
for evolving the spacetime metric is substantively the 
sa me as the generalized harmonic method described 
by iPretoriusI ()2005bl lal). except that the FD scheme is 
fourth order accurate and uses fourth order Runge- 
Kutta time integration. We model the NS material 
as a perfect fiuid and solve the hydrodynamic energy- 
momentum and continuity equations using conserva- 
tive high-resolution shock-capturing schemes with sec- 
ond order Runge-Kutta time integration, and enforce 
strict conservation at AMR boundaries using the fiux 
correction method of iBerger Golellal (|T989h . We 
have implemented severa l metho ds for calculating inter- 
cell fiux es fHLL. nHarten et al.l [1983: the Roe solver, 
[Eulderin k fc Mellemal 119951 : and the Marquina fiux, 
iMarquina et al.lll992f ) and for reconstructing fiuid prim- 
itive variabl es at cell interfaces (MC and minmod. iTorol 
T997[ PPM. IGole lla fc Woodward! [1981 and WENO-5, 
Tchekhovskoy et al...2007.) . Unless otherwise noted, the 



simulations described below were performed with HLL 
and WENO-5. Our hydrodynamics scheme allows for 
any EOS of the form P = P(p, e) (e.g., F-law, piecewise 
polytrope, and tabular EOSs). We have tested the new 
hydrodynamics sector of our code on problems includ- 
ing ID and 2D Riemann problems, Bondi accretion, and 
single NSs. More details on our c ode and tests will be 
presented in [Stephens "eFall (|20TT[ ). 

For this first study the only parameter we vary is the 
initial periapsis separation. The BH and NS have a 4:1 
mass ratio, and both are initially non-rotating. The NS 
is modeled as a TOV star with a broken F-law EOS (la- 
beled "HB" in lR.ead et al.l 120091 and including a ther- 
mal component to allow for shock heating) and has mass 
1.35M0 and radius 11.6 km. The initial orbital param- 
eters describe a hyperbolic encounter with relative ve- 
locity w = 1000 km s~^, correspon ding to the central 
region of a NC (jO'Learv et al.ll2009[ ). These orbits are 
nearly parabolic, with e — 1 ^ O(10~^), and hence the 
close-encounter behavior also adequately describes such 
events in GCs. We superimpose initial data for the BH 
and NS at an initial separation of 50M (498 km) with 
initial velocities according to a Newtonian orbit with the 
desired Vp. Though these superposed initial data do not 
strictly satisfy the constraint equations except at infinite 
separation, tests performed at various initial separations 
indicate that 50M is sufficiently large that the constraint 
violation does not appreciably affect the system (relative 
to truncation error). 

3. RESULTS AND DISCUSSION 

We consider a range of periapsis separations from 
Tp/M = 5.0 to 15 {i.e., 50 to 150 km). (Hence- 
forth, we will consider Vp to be normalized by M.) 
In all of these cases, sufficient energy is carried away 
by GWs to result in a bound system. Our simu- 
lations exhibit three types of behavior: (1) a direct 
plunge {rp = 5.0,5.83,6.67,6.81), (2) following the ini- 
tial periapsis passage, a single elliptical orbit and then 
a plunge {vp = 6.95,7.22,7.5), and (3) following the 
initial periapsis passage, a long-period elliptical orbit 
{rp = 8.75,10.0,12.5,15.0). For the latter group (and 
the high resolution rp = 7.5 run), the entire orbit is pro- 
hibitively long to simulate, and we focus on the burst of 
GWs associated with the first periapsis passage. For one 
case in each class {rp = 5.0, 7.5, 10.0) we ran three sim- 
ulations with different characteristic mesh spacings (but 
always with 7 refinement levels) for convergence studies. 
At t = 0, the low (medium, high) resolution run had 
finest meshes covering the BH and NS of roughly 80^ 
(100^,150^) cells, resolving the NS diameter with ~ 40 
(50, 75) cells and the BH horizon diameter with ~ 70 
(85, 130) cells. (We note that the level structure is set 
by truncation error estimates and is adjusted with time.) 

All other simulations were run at medium resolu- 
tion. Unless otherwise noted, results will be reported 
for medium resolution, with error bars (where appropri- 
ate) computed from convergence calculations. Our sim- 
ulations employ compactified coordinates such that the 
outer boundaries extend to spatial infinity. Thus, the 
global (ADM) M and J should be conserved. In prac- 
tice, however, we must evaluate these quantities at a fi- 
nite distance, making them subject to gauge artifacts, 
some propagating outward from the central BH/NS re- 
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gion from t = 0. For t < 200M, an extraction sphere 
of 300M is free of propagating artifacts, whence M (J) 
is conserved to better than 0.3 (2.0)% for ah cases at 
medium resolution. 

Based on the abov e of runs, som e guidance from per- 
turbative results (Peters fc Mathews 1963 : Turner, J 9771 : 
iBerrv Ga'iill2010D . and a zoom- whirl geodesic analogu e 
of the two body problem (jPretorius fc Khuranal [2QQ7I ) . 
we conjecture the following qualitative behavior as a 
function of Tp for initial encounters resulting in a bound 
system. 

Consider n, the no n- negative, integer number of peri- 
apsis passages before disruption/plunge, and , the peri- 
apsis distance on the i^^ encounter for 1 < i < n. Define 
r^+-^ to be an effective periapsis distance for the final 
close encounter (i.e., the corresponding before the dis- 
ruption/plunge part of the final encounter). Group (1) 
above has n = 0, group (2) n = 1, and group (3) n > 1. 
The behavior of a close encounter will depend sensitively 
on the distance Sr^ = Tp—Vc between and a radius Vc of 
an effective unstable circular orbit with the same energy 
and angular momentum. If ^r* is sufficiently small (rela- 
tive to M), the orbit will exhibit a whirl phase, where it 
asymptotes to a nearly circular orbit. The smaller ^r* , 
the longer the duration of the whirl, with a maximum 
when 6rp = that equals the time required for the bi- 
nary to lose its excess orbital kinetic energy, either via 
GW emission or tidal transfer of energy to the NS ma- 
terial. The excess energy is the difference between the 
total energy of the binary entering the whirl phase and 
a putative binary on a quasi- circular inspiral at Tc. 

Because of the requirement that (5r* be small, we only 
expect the possibility of significant whirling near the ul- 
timate or penultimate encounter. If Sr^ is negative, the 
whirl will directly transition to a plunge. If positive and 
small, there will be a separation following the whirl; how- 
ever the effective Tc for the next encounter increases while 
decreases due to GW emission, and since this is quite 

sizeable for our 4:1 mass ratio system, (5r*+^ will likely be 
negative, resulting in a subsequent plunge. Furthermore, 
the separation where the NS starts to be tidally disrupted 
is within the radii of unstable whirl orbits, hence the 
NS will not survive any prolonged whirl phase. Prior 
encounters (for larger n cases) will not exhibit signifi- 
cant whirling, and while (5r* is large the orbital evolution 
could better be described as a series of precessing ellipses 
with decreasing eccentricity and semi- major axis. 

Considering the number of orbits n as a function of the 
initial rp^ n{rp) is monotonically increasing, and the val- 
ues of rp where one could see a notable final whirl would 
be near the steps in n{rp). The most pronounced whirl 
behavior will occur for small n{rp)] as n{rp) increases, 
the amount of excess kinetic energy left over once the 
final orbit is reached decreases, and for sufficiently large 
n the late stages will essentially follow a quasi-circular 
inspiral. 

Figure [T] illustrates some of the varied phenomena en- 
countered near the transition between n{rp) = and 
n{rp) = 1, occurring at rp ^ 6.88 ± 0.08. Most striking 
is the amount of rest mass remaining after the merger 
as a function of rp (Table [T]). For cases in which the NS 
plunges without significant disruption, such as = 5 or 



rp = 7.5, less than 1% of the initial mass is available to 
form an accretion disk. However, for rp closer to the first 
transition (r^ = 6.67 and 6.81, see Figure [T]), the NS is 
stretched into a long tidal tail, and a sizeable amount of 
bound material is left to form an accretion disk — 12% of 
the initial NS rest mass for rp = 6.81. 

Figure [2] shows the approximate rate of fallback as 
a function of time for rp = 6.81, 6.95, 7.22, and 7.5. 
This is the rate at which material on elliptical orbits is 
expec ted to return to the accretion disk (see iRosswogI 
120071 ). (These accretion rates are likely upper lim- 
its since they do n ot account for nuclear burning, see 
iMetzger et al.|[2QTQl ) The fallback rate is larger for cases 
with larger disk masses, such as rp = 6.81, but all cases 
exhibit an approximate t~^^^ falloff. This time depen- 
dence was predic ted fo r stella r disruptions around super- 
massive BHs bv lRi^ (dill). It appears unlikely that 
BH-NS mergers with fallback rates as in Fig. [2] will be 
able to explain sGRBs with extended emission (see, e.g., 
iNorris fc Bonnell l2006l ) if this emission is due to feed- 
ing of the accretion disk at late times. For example, by 
t ^ 100 s, the luminosity for the = 6.81 case would be 

only L ~ rjMc^ ~ 2 x 10^^ erg/s (assuming an efficiency 
7/ = 0.1). 

Table [T] also shows the total energy and angu- 
lar momentum lost to GWs for rp between 5 and 
12.5. For the cases that we followed through merger, 
we find 0.7-1.7% of the total mass lost to GW en- 
ergy, and estimate the final spins of the BHs to 
be (0.49 ± 0.01,0.45,0.37,0.47,0.50,0.50) for rp = 
(5.00, 6.67, 6.81, 6.95, 7.22, 7.50) respectively. The energy 
loss is largest for the transitional case (r^ = 6.95), which 
has a large pulse from the whirling first passage and a 
second burst from the merger (Fig. [3|). Table [T] also 
shows the GW losses for the initial encounter in cases 
where the NS survives the periapsis passage (columns 
5 and 6). These fiy-by p ulses can be compared with 
the prediction of iTurnerl (|1977l ) (Fig. [3l lower panels, 
and Fig. [4]), who used the Newtonian orbit together with 
quadrupole physics for the GW emission, which we will 
call the Newtonian Quadrupole Approximation (NQA). 
Our waveforms show roughly the same pulse shape as 
the NQA prediction but have larger amplitudes for the 
smaller rp cases. At = 15 we find the gravitational 
waveform from the initial fiy-by to be indistinguishable 
from the NQA prediction at our resolution (±10%). 

The enhancement in GW energy losses for close en- 
counters may be due (in part) to zoom-whirl-like behav- 
ior. Figure [1] shows the GW energy loss as a function of 
Tp, along with the NQA pr ediction and a fit consistent 
with zoom- whirl dynamics (jPretorius fc Khuranal [2QQ7I ) 
in the regime (r^ < 10) where we start to see significant 
departures from the NQA approximation. 

4. CONCLUSIONS 

An interesting result of this general relativistic study 
that is quali tatively co nsistent with previous Newtonian 
studies (e.g. iLee et af] 12010). is the great variability of 
the outcome as a function of impact parameter. For ex- 
ample, the remnant disk masses following merger range 
from nearly zero up to ~ O.3M0. Mergers leading to sig- 
nificant disks occur in a small (but not negligible) region 
of parameter space and could produce a sGRB. 
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Fig. 1. — Rest mass density in the equatorial plane from BH-NS simulations with varying rp. The four panels on the left show (left to 
right, top to bottom, same color scale): (1) the BH and NS merging (t = 6.82 ms, = 5), (2) the NS being stretched into a long tidal 
tail (t = 9.97 ms, rp = 6.81), (3) a brief mass transfer episode during the NS's first periapsis passage (t = 8.78 ms, rp = 6.95), and (4) the 
NS's subsequent distortion (t = 10.6 ms, rp = 6.95). On the right is a nascent accretion disk (top, t = 12.3 ms, rp = 6.67), and a late-stage 
accretion disk (bottom, t = 38.2 ms, rp = 7.5, low resolution, PPM). The color scale is logarithmic with units of the initial maximum 
density (po = 8.3 X 10-*^^ g cm"^). The black hole is roughly the same coordinate size in all panels {Rbu = 16 km), which can be used to 
infer the relative scale of each snapshot. 
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TABLE 1 

Disk properties and GW energy and angular momentum losses 



^Rest mass remaining outside the BH shortly after merger, normalized by the initial total rest mass. 
^Unbound rest mass estimated using local fluid velocities. 

^Rough initial accretion timescale (race = Mq/Mq) evaluated shortly after merger. 

•^Energy and angular momentum lost to GWs during the flrst close encounter. 

^ Total GW energy and angular momentum losses for cases which were followed through merger. 

^Results are from medium resolution runs; values in parentheses are Richardson extrapolated estimates using low and high resolutions, 
where available. Note that the relatively large error for rp = 7.5 (and to a lesser extent rp = 5, 10) is due in part to truncation error 
altering the actual periapsis by a small amount, and in this regime the GW emission is highly sensitive to binary separation (Figure 0J. 



In follow up work we plan to extend the parameter 
space survey, varying the NS EOS and BH spin. We 
also intend to explore the detectability of these events 
with GW detectors, and conclude with brief prelimi- 
nary comments. These signals may be difficult to de- 
tect with instruments such as LIGO since they lack a 
long inspiral phase and most of the power is at high fre- 
quencies (1500-2000 Hz for the masses considered here). 
Using the broadband AdLIGO noise curve, we find sky- 
averaged SNR of 3-9 for = 5 — 10, assuming a dis- 
tance of 100 Mpc. Scaling the GW emission up to 
(Mns,^bh) = (2,8)M0 and assuming optimal orien- 
tation gives SNR of 8 out to 340 Mpc for Tp = 7.5. En- 



counters with larger Tp could be easier to detect since 
they will have a number of fly-by GW bursts before 
the merger. Given that some BH-NS capture binaries 
emit the majority of their GW power at the high fre- 
quency end of the LIGO noise curve and could be rela- 
tively weak (both in terms of GW emission and extrap- 
olated electromagnetic emission based on disk mass), it 
may be a wor thwhile exercise to revisit the analysis of 
GRB 070201 (jAbbott et al.ll2008[ ) with burst-templates 
adapted to capture driven BH-NS encounters. 

We thank Adam Burrows, John Friedmann, Roman 
Gold, Benjamin Lackey, and Richard O'Shaughnessy 
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Fig. 2. — Approximate fallback accretion rates for rp = 6.81, 
6.95, 7.22, and 7.5. These rates are evaluated at times ranging from 
1.0 to 2.2 ms after the approximate time when the BH accretion 
rate plateaus following the merger. For this diagnostic, we consider 
the fluid in each cell as a ballistic particle and take its orbital 
period as the approximate fallback timescale. The instantaneous 
BH accretion rates evaluated at the same time are shown at the 
upper right (arbitrary abscissa). 
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-1 1 

time (ms) 

Fig. 3. — (upper panel) The Newman-Penrose scalar ^4 on the z- 
axis (orthogonal to the orbit) for rp = 6.95. The flrst pulse is from 
the initial close encounter, the second from the merger-ringdown. 
Between the pulses there is an oscillation due to the rotating, dis- 
torted neutron star, which is signiflcantly torqued during the flrst 
encounter. Here t = corresponds to the start of the simula- 
tion, (lower panels) The real and imaginary components (black 
diamonds and red squares) of the I = 2, m = 2 spherical harmonic 
of r^4 for rp = 7.5 (left) and rp = 10 (right). For comparison the 
NQA analytical results are shown multiplied by an overall factor 
so that the magnitude and phase match at peak (t = 0). 
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